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Abstract:  Eoam  fraction  can  be  retrieved  from  space-based  microwave  radiometric  data  at 
frequencies  from  1  to  37  GHz.  The  retrievals  require  modeling  of  ocean  surface  emissivity 
fully  covered  with  sea  foam.  To  model  foam  emissivity  well,  knowledge  of  foam 
properties,  both  mechanical  and  dielectric,  is  necessary  because  these  control  the  radiative 
processes  in  foam.  We  present  a  physical  description  of  foam  dielectric  properties  obtained 
from  the  foam  dielectric  constant  including  foam  skin  depth;  foam  impedance;  wavelength 
variations  in  foam  thickness,  roughness  of  foam  layer  interfaces  with  air  and  seawater;  and 
foam  scattering  parameters  such  as  size  parameter,  and  refraction  index.  Using  these,  we 
analyze  the  scattering,  absorption,  reflection  and  transmission  in  foam  and  gain  insights 
into  why  volume  scattering  in  foam  is  weak;  why  the  main  absorption  losses  are  confined  to 
the  wet  portion  of  the  foam;  how  the  foam  impedance  matching  provides  the  transmission  of 
electromagnetic  radiation  in  foam  and  maximizes  the  absorption;  and  what  is  the  potential 
for  surface  scattering  at  the  foam  layers  boundaries.  We  put  all  these  elements  together  and 
offer  a  conceptual  understanding  for  the  high,  black-body-like  emissivity  of  foam  floating  on 
the  sea  surface.  We  also  consider  possible  scattering  regimes  in  foam. 
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1.  Introduction 

Oceanic  whitecaps  are  a  conspicuous  expression  of  wave  breaking  with  air  entrainment.  Creating 
rafts  of  sea  foam  on  and  bubble  plumes  below  the  sea  surface,  breaking  waves  transfer  momentum, 
heat,  and  mass  across  the  air-sea  interface  [1-4]  and  alter  the  optical  field  of  the  water  column,  the 
surface  roughness,  and  the  ambient  noise  in  the  ocean  [5-7].  One  (among  others)  physical  quantity  that 
is  used  to  evaluate  these  air-sea  interaction  processes  is  the  whitecap  fraction,  defined  as  the  fraction  of 
the  ocean  surface  covered  by  foam. 

Whitecap  fraction  can  be  evaluated  with  a  number  of  parameterizations  developed  for  open 
ocean  [8]  and  coastal  zone  [9]  using  photographic  measurements  [10-12]  in  the  visible  portion  of  the 
electromagnetic  (EM)  spectrum.  Whitecaps  can  also  be  observed  at  infrared  [13]  and  microwave 
wavelengths  [14].  Recent  efforts  to  retrieve  whitecap  fraction  from  space-based  radiometric 
observations  of  the  ocean  surface  at  microwave  frequencies  from  1  to  37  GHz  show  promise  [15-18]. 
The  retrieval  algorithm  requires  the  modeling  of  the  ocean  surface  emissivity  fully  covered  with  sea 
foam  ef. 

Laboratory  and  field  measurements  have  shown  that  e/  at  microwave  frequencies  is  high,  close  to 
that  of  a  blackbody  [14,19-22].  To  explain  and  model  e/,  one  needs  to  understand  the  radiative 
processes  taking  place  in  foam.  Knowledge  of  the  mechanical  and  dielectric  properties  of  sea  foam  is 
necessary  for  such  an  understanding.  While  extensive  oceanographic  research  facilitates  the 
characterization  of  the  foam  mechanical  properties  (Chapters  4.4  and  4.5  in  [23]),  the  knowledge  of  the 
foam  dielectric  properties  is  incomplete. 

Basis  for  the  dielectric  properties  of  a  medium  is  its  relative  dielectric  constant  (permittivity). 
Permittivity  determines  such  intrinsic  characteristics  of  the  medium  as  skin  (and  penetration)  depth, 
impedance,  refractive  index,  etc.  These  quantities  control  the  attenuation  due  to  absorption  and  scattering 
and  the  transmission  and  reflection  of  the  EM  radiation  in  the  medium.  These  radiative  processes  affect, 
in  turn,  the  energy  lost  in  and  then  emitted  by  the  medium.  Therefore,  an  understanding  of  the  relative 
importance  of  all  these  characteristics  and  processes  can  assist  the  choice  of  an  approach  to  model  e/,  the 
formulation  of  relevant  assumptions  and/or  simplifications,  and  the  interpretation  of  the  model  results. 

We  address  the  need  for  more  complete  characterization  of  foam  dielectric  and  radiative  properties 
in  a  series  of  papers.  Anguelova  [24]  analyzed  the  available  information  to  determine  suitable  formula 
to  predict  the  complex  permittivity  of  sea  foam  e/.  Anguelova  and  Gaiser  [25]  focused  on  the  skin 
depth  of  foam  layers  with  a  vertical  void  fraction  profile.  In  the  present  paper,  we  consider  the  foam 
properties  at  microwave  frequencies,  which  control  fundamental  radiative  processes  such  as  reflection, 
scattering,  and  transmission  in  vertically  structured  foam  layers.  The  parameters  used  for  this 
consideration  are  foam  impedance  rj,  size  parameter  x,  and  refractive  index  m  (Section  3).  Analysis  of 
these  dielectric  properties  helps  to  identify  the  unique  traits,  which  the  fundamental  radiative  processes 
acquire  when  microwave  radiation  interacts  with  sea  foam  (Section  4).  Einally,  we  put  all  our  findings 
together  to  formulate  a  concept  of  the  high  foam  emissivity  (Sections  5.1  and  5.2)  and  discuss 
scattering  regimes  in  foam  (Section  5.3).  The  information,  findings,  and  generalizations  in  this  series 
of  papers  form  the  physical  basis  on  which  we  have  built  our  foam  emissivity  model  (briefly  described 
in  [16]  and  to  be  detailed  in  a  forthcoming  paper). 
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2,  Background 

2.1.  Sea  Foam  as  a  Medium 

Both  whitecaps  on  the  surface  and  bubble  plumes  beneath  them  constitute  sea  foam.  For  passive 
remote  sensing  of  the  ocean  surface,  the  foam  skin  depth  dictates  our  interest  in  the  floating  foam 
layers,  excluding  deeper  bubble  plumes  [25].  Whitecaps  corresponding  to  the  inception  of  wave 
breaking  are  usually  thick  and  bright  (visually  and  radiometrically)  and  are  referred  to  as  active 
whitecaps.  Whitecaps  in  their  decaying  phase  are  thinner  and  dimmer  and  are  referred  to  as  residual 
whitecaps.  Anguelova  [24]  gives  an  extended  definition  of  the  sea  foam. 

Sea  foam  in  surface  layers  has  a  specific  vertical  structure  comprising  large  thin-walled  bubbles 
close  to  the  air-foam  interface  characterized  with  high  air  content  (dry  foam)  and  smaller,  thick-walled 
bubbles  with  high  water  content  (wet  foam)  close  to  the  seawater  boundary  [20].  A  set  of  macroscopic 
quantities  such  as  void  fraction/,  (defined  as  the  fraction  of  a  unit  volume  of  seawater  that  is  occupied 
by  air)  and  foam  layer  thickness  t  can  describe  this  mechanical  structure.  A  set  of  microscopic 
quantities  such  as  bubble  dimensions  (radius  a  and  wall  thickness  w)  and  concentration  or  size 
distribution  N{a)  are  useful  when  one  considers  the  characteristics  of  the  bubbles  forming  the  foam. 

To  describe  oceanic  whitecaps  under  various  conditions  and  in  different  lifetime  stages,  the  full 
range  of  void  fractions  fa  (from  0  to  1)  needs  to  be  considered  [24].  It  was  shown  that  various 
functional  forms  could  represent  the  shape  of  the  void  fraction  profile  in  the  foam  depth  [25].  A  review 
of  oceanographic  data  from  laboratory  and  field  experiments  [25]  established  a  plausible  range  of  foam 
layer  thicknesses  encountered  in  the  ocean,  from  1  cm  to  more  than  12  cm  in  active  whitecaps  and 
from  0.1  to  1  cm  when  the  whitecaps  decay. 

2.2.  Bubbles  in  a  Foam  Layer 

Bubble  size  distributions  that  characterize  bubble  plumes  at  different  depths  and  under  various  wind 
speeds  have  been  published  extensively;  see  summaries  in  [23,26,27].  Meanwhile,  there  are  just  a  few 
reports  for  bubble  size  distributions  within  the  foam  layer.  Militskii  et  al.  (Figure  3  in  [28])  report 
bubble  diameters  within  the  foam  layer  in  the  range  from  0.2  mm  to  2  mm  with  the  most  probable 
value  being  0.3  mm  to  0.5  mm.  From  analysis  of  20  foam  images,  Guo  et  al.  [29]  report  the  mode  of 
the  bubble  size  distribution  to  be  around  0.5  mm  “in  the  bottom  half  of  the  foam  layer”.  Certainly,  the 
range  of  bubble  sizes  is  much  wider  in  natural  conditions  under  various  wind  speeds.  To  encompass  as 
fully  as  possible  this  range  of  bubble  sizes,  we  refer  to  recent  laboratory  and  field  measurements  of 
bubble  size  distributions  in  bubble  plumes  close  to  the  surface,  about  0.5  m  below  the  interface,  and 
soon  after  the  breaking  event  [27,30].  Thus,  we  consider  bubbles  with  radii  a  of  the  air  cavity  from 
0.05  mm  to  10  mm  with  a  broad  peak  around  radii  of  0.05  mm  to  1  mm. 

2.2.  Definitions  of  Electromagnetic  Properties 

Basic  physical  relationships  from  the  electromagnetic  theory  determine  the  dielectric  properties  of  a 
medium.  Table  1  lists  the  formulae  used  to  obtain  all  foam  dielectric  properties  from  the  complex 
permittivity  of  foam  £/at  microwave  frequencies  F  from  1.4  to  37  GHz.  The  formulae  used  are  valid 
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for  scatter  free  medium,  but  are  applieable  to  sea  foam  beeause  a  survey  of  experimental  and  modeling 
observations  has  shown  that  at  the  frequeneies  of  interest  the  seattering  in  foam  is  weak  to  negligible 
(Seetion  2.3  in  [24]). 

Diseernible  interfaees  at  the  air-foam  and  foam-water  boundaries  point  toward  possible  eontribution 
to  radiation  attenuation  due  to  surfaee  seattering  by  foam.  To  investigate  whether  the  foam  interfaees 
with  air  and  water  are  rough  or  smooth  as  seen  by  mierowave  radiation,  we  use  two  eriteria  appropriate 
at  mierowave  frequeneies  ([32],  p.  827).  The  Fraunhofer  eriterion  defines  a  surfaee  as  smooth  if 
a  <  Oref,  where  (7  is  the  standard  deviation  of  the  surfaee  height  (or  RMS  height)  and  the  referenee 


32^- 


with  A,  being  the  wavelength  and  6  the  ineidenee  angle  of  the  EM  radiation. 


The  other  eriterion  expresses  <7  in  eleetromagnetie  units  with  the  quantity  k<7  where  k  =  2tiIX  is  the 
radiation  wavenumber. 


3,  Foam  Dielectric  Properties 

3.1.  Sea  Foam  Permittivity 

Systematic  investigation  of  various  mixing  rules  yielded  several  possibilities  for  computing  the 
dieleetrie  constant  of  sea  foam  £/at  microwave  frequencies  [24].  For  consistency  with  the  study  on  the 
foam  skin  depth  [25],  here  we  use  £/obtained  with  the  Polder-van  Santen  (PS)  mixing  rule  (Equation  (1) 
in  Table  1).  The  calculations  use  a  double  Debye  model  for  the  complex  dielectric  constant  of  seawater 
£•[33].  Refer  to  [24]  for  the  dependencies  of  obtained  with  the  PS  mixing  rule  on  the  frequency,  sea 
surface  temperature  (SST),  and  salinity.  All  results  are  illustrated  with  an  exponential  profile  for  the 
foam  void  fraction  y^(z)  with  values  at  the  air-foam  and  foam-water  boundaries  of  99%  and  1%, 
respectively.  Refer  to  [25]  for  comparison  of  exponential  fa{z)  to  other  functional  forms  of  the  void 
fraction  profile. 

3.2.  Foam  Skin  and  Penetration  Depths 

Anguelova  and  Gaiser  [25]  obtained  and  analyzed  in  detail  the  skin  depth  d  (Equation  (2)  in  Table  1) 
of  sea  foam  vertically  stratified  in  layers  with  thickness  t  from  0.2  cm  to  10  cm.  For  various 
frequencies  F  and  the  chosen /,(z)  profile  (Section  3.1),  d  could  vary  from  0.17  cm  to  no  more  than 
7  cm.  Analysis  of  the  relationship  between  d,  t,  and  the  fa  range  (upper  to  lower  values  at  the  foam 
layer  boundaries)  showed  that  the  thermal  emission  from  foam-covered  surfaces  could  be  from  the 
entire  foam  layer,  from  part  of  the  foam  layer,  or  from  both  the  foam  layer  and  the  seawater. 

For  any  frequency  or  any  layer  thickness,  d  is  thicker  than  its  corresponding  penetration  depth  S 
(Equation  (3)  in  Table  1)  by  a  factor  of  up  to  2  [25].  Figure  1  shows  tJas  a  function  of  frequency  for 
three  fixed f,  values.  For/,  =  90%,  S ranges  from  a  few  centimeters  to  a  few  meters;  for f,  <  60%,  S is 
mostly  less  than  1  cm. 
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Table  1.  Basic  relationships  from  the  eleetromagnetic  theory  used  to  obtain  foam  dieleetrie  properties. 

Equation  # _ Property _ Symbol  (Units) _ Formula _ Reference _ Notes _ 

Dielectric  properties 

Seawater  is  environment  with 
permittivity  e. 

Bubbles  are  inclusions  with 
void  fraction  fa 

d 


cc{z)dz  —  1  a  field  attenuation  coefficient 


(2) 

Skin  depth 

d  (mm) 

0 

a{z)-  ^  '  .Im^£-^(z) 

p.  847  in  [32] 

F  frequency  (Hz) 
c  speed  of  light  (cm  s  ') 

(3) 

Penetration  depth 

^(mm) 

s 

1  k^{z)dz  =  1 

0 

d  =  dl2  @  e^  {z)  =  const 

ditto 

Scattering  ignored 
Extinction  =  Absorption 
i.e.,  k^=k^=  2a 

(4) 

Intrinsic  impedance 

V 

77  =  1/^ 

p.  226  in  [35] 

Normalized  (relative),  complex 

(5) 

Wavelength  in  foam 

Xf- (cm) 

II 

p.  1453  in  [36] 

Xq  free-space  wavelength 

(6) 

Propagation 
constant 
(wave  number) 

kf{cvn^) 

kj.  =  iTijXf 

Scattering  parameters 

p.  116  in  [37] 

(7) 

Size  parameter 

X 

x  =  a  -  ylnj j  =  k^a 

p.  128  in  [37] 

a  bubble  radius 

(8) 

Refraction  index 

m 

II 

II 

p.  116  in  [37] 

m'  and  m"  real  and  imaginary 
parts  of  m 

Dielectric  constant 


£  2^*  +  2(£ y 


-  f  1-g 

l  +  le  +  lie, 


Equation  (9.7)  in  [34] 
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Figure  1.  Foam  penetration  depth  as  a  funetion  of  frequeney  d  {F)  obtained  with 
Polder-van  Santen  mixing  rule  with  exponential  void  fraetion  profile  ranging  from  99%  at 
the  air- foam  interfaee  and  1%  at  the  foam-seawater  boundary  at  fixed  seawater  temperature 
{Ts  =  20  °C)  and  salinity  {S  =  34  psu)  at  three  values  of  foam  void  fraetion  fa. 


ioy 


3.3.  Foam  Impedance 

The  eoneept  of  impedanee  has  been  applied  in  different  fields  of  physios  ([38];  p.  282).  In  the  EM 
theory,  the  intrinsio  impedanee  of  a  medium  follows  from  the  relationship  between  the  amplitudes  of 
the  eleotrio  and  magnetio  fields  of  EM  radiation.  It  depends  only  on  the  medium  properties  and  is 
important  when  the  propagation  of  EM  radiation  through  a  medium  is  oonsidered. 

Figure  2.  Normalized  intrinsio  impedanee  of  sea  foam  obtained  from  foam  permitivity 
(Equation  (4)  in  Table  1)  at  fixed  seawater  temperature  {Ts  =  20  °C)  and  salinity  {S  =  34  psu): 

(a)  as  a  funetion  of  frequeney  for  seawater,  air,  wet  foam  {fa  =  10%),  and  dry  foam 
{fa  =  98%);  (b)  as  a  funetion  of  foam  void  fraetion  at  three  frequenoies. 


Eigure  2(a)  shows  the  magnitude  of  a  oomplex  foam  impedanee  rj  (Equation  (4)  in 

Table  1)  for  dry  (oonstant  vertioal  profile  of  /,  =  98%,  dash-dot  line)  and  wet  (eonstant  profile  of 
fa  =  10%,  dashed  line)  foam.  The  impedanees  of  the  seawater  (solid  line)  and  air  (dotted  line)  are 
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shown  for  comparison.  Figure  2(b)  depicts  how  the  foam  impedance  changes  as  a  function  of/,  in  a 
vertically  structured  foam  layer  for  three  frequencies.  Over  the  range  of  considered  frequencies,  foam 
impedance  varies  from  1  to  less  than  0.2  as/,  decreases  from  100%  to  0%. 

3.4.  Wavelength  Changes  in  Foam  Layers 

Vertical  variations  of  the  real  part  of  the  permittivity  in  foam  thickness  invoke  changes  in  the 
radiation  wavelength  //and  propagation  constant  k/as  compared  to  that  in  air  Aq  (Equations  (5)  and  (6) 
in  Table  1).  Table  2  shows  values  of  //in  dry  and  wet  foam  represented  with  /  =  98%  and  10%, 
respectively,  for  all  considered  frequencies.  The  corresponding  values  in  seawater  A,  are  also  included 
for  reference. 

Table  2.  Wavelengths  at  fixed  SST  (20  °C)  and  salinity  (34  psu)  in  air  (//),  in  seawater 
(A),  in  dry  and  wet  foam  (//). 


FGHz 

Ao  cm  (in  air) 

Af 

cm 

A  cm  (in  seawater) 

fa  =  98% 

fa=10Vo 

1.4 

21.4 

20.8 

2.762 

2.55 

6.8 

4.4 

4.28 

0.603 

0.56 

10.7 

2.8 

2.72 

0.414 

0.38 

18.7 

1.6 

1.56 

0.286 

0.26 

23.8 

1.3 

1.22 

0.253 

0.23 

37.0 

0.8 

0.79 

0.207 

0.19 

Figure  3.  Wavelength  in  foam  //(cm)  (Equation  (5)  in  Table  1)  as  a  function  of  foam  void 
fraction  at  fixed  seawater  temperature  (/  =  20  °C)  and  salinity  (5*  =  34  psu)  and  different 
frequencies.  Vertical  lines  divide  the  /,  range  into  three  regions  whose  values  can  be 
associated  with  bubble  diamaters  predominant  for  each  region. 


Eigure  3  depicts  the  change  of  //  from  //  to  /I  as  it  propagates  through  foam  with  varying  void 
fraction.  In  the  figure,  for  each  frequency  we  discern  three  regions  in  the  trend  of  //based  on  the  rate 
of  change  dXJdfa.  In  dry  foam,  //varies  slowly  remaining  close  to  //;  dXJdfa  changes  by  a  factor  of 
1.4  (1.1)  at  1.4  (37)  GHz  as  void  fraction  changes  from  80%  to  100%.  Eor  void  fractions  below  40%, 


Remote  Sens.  2012,  4 


1169 


Afis  quite  close  to  A,  and  also  changes  relatively  slowly.  The  slope  dXjIdfa  changes  by  a  factor  of  2.6 
on  average  over  the  range  of  frequencies  when  varies  from  40%  to  10%.  A  steeper  transition 
between  the  air-like  and  seawater-like  behavior  takes  place  in  the  void  fraction  range  of  80%  to  40% 
with  d  XjJd  fa  changing  on  average  by  a  factor  of  4.8.  The  lower  the  frequency  the  more  clearly  these 
three  regions  are  seen. 

On  this  basis  we  can  broadly  represent  a  vertically  structured  foam  layer  as  a  sequence  of  three 
relatively  homogeneous  sublayers,  which  stand  for  dry,  intermediate  wet,  and  wet  foams.  We  identify 
these  sublayers  by  their/,  values  with  vertical  lines  in  Figure  3.  While  we  base  this  conceptual 
representation  on  different  rates  of  change  of  /^values,  Raizer  (Figure  1(a)  in  [39])  gives  a  similar 
schematic  of  a  vertically  stratified  foam  based  on  experimental  observations  [28]. 

3.5.  Scattering  Parameters  of  Foam 

3.5.1.  Characteristic  Bubbles  for  Each  Foam  Sublayer 

The  known  scattering  rules  of  individual  particles  [37]  can  be  applied  to  foam  mixture  because  the 
foam  macroscopic  parameters  (e.g.,/,  and  Af  are  related  to  those  of  a  separate  bubble  {e.g.,fb  and  /I*). 

Qualitatively,  such  a  relationship  is  possible  because  the  vertical  structure  of  a  foam  layer  is 
associated  with  segregation  by  size  of  the  bubbles  forming  the  foam  (see  [24]  and  Section  2.1).  Indeed, 
we  expect  that  dry  foam  consist  predominantly  of  large  bubbles  with  an  upper  limiting  radius  of 
10  mm  (Section  2.2).  Meanwhile,  most  of  the  smaller  bubbles  with  a  lower  limiting  radius  of  0.05  mm 
would  concentrate  preferentially  in  the  lower,  wet  part  of  the  foam  layer.  The  intermediately  wet  foam 
between  these  two  limiting  portions  would  comprise  bubbles  with  wider  size  variation,  with  most 
typical  sizes  around  the  peak  of  the  bubble  size  distribution  (Section  2.2).  Thus,  the  three  sublayers  of 
vertically  structured  foam  associated  with  specific  /,  values  (Section  3.4)  could  be  alternatively 
differentiated  by  expected  specific,  predominant  bubble  radius  a  (Figure  3). 

Quantitatively,  this  relationship  is  based  on  the  fact  that  foam  void  fraction  can  be  obtained  from 
the  bubble  size  distribution  [7].  The  void  fraction  /,  of  a  collection  of  A  bubbles  with  a  fixed  outer 
radius  a  and  a  wall  (shell)  of  seawater  w  =  a-qa,  where  qa  is  the  inner  radius  of  a  bubble,  is 

f  ^  f  _  1  v  air  volume  in  a  bubble  _  3 
"  ~  *  A  “  total  volume  of  a  bubble  ^ 

For  example,  if  a  =  10  mm  and  w  =  100  pm  {i.e.,  q  =  I  -  wla  =  0.99),  then  f  ~  97%.  Since  a  bubble 
with  such  dimensions  (large  a,  thin  wall  w)  is  characteristic  for  the  dry  foam  sublayer  represented  with 
fa  =  98%,/,  can  characterize  the  mixture’s  fa.  Similarly,  the  void  fractions  of  smaller  bubbles,  e.g., 
a  =  0.3  mm  and  0.1  mm,  having  walls  with  a  thickness  w  of  0.05  mm  (i.e.,  q  =  0.83  and  0.5  respectively) 
are  representative  of  foam  with/,  =fb  =  q  ~  60%  and  12%,  respectively. 

3.5.2.  Size  Parameter  in  a  Foam  Layer 

The  size  parameter  x  (or  diffraction  parameter  [40])  describes  the  relationship  between  the  radiation 
wavelength  in  a  mixture  and  the  sizes  of  the  scatterers  in  it.  In  the  case  of  foam,  x  quantifies  how  this 
relation  changes  in  foam  layer  thickness  as  both  Af  and  the  dominant  bubble  size  a  change. 
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Figure  4  shows  a  family  of  curves  for  all  considered  frequencies  representing  x  as  a  function  of^  of  a 
foam  mixture  comprising  bubbles  with  a  fixed  radius  a  =  0.3  mm  (a  size  within  the  peak  of  the  bubble 
size  distribution,  Section  2.2).  Changes  of^^  in  the  depth  of  foam  comprising  bubbles  with  the  fixed 
radius  a  are  realized  via  variations  of  bubble  wall  thickness  w  which  yield  q  =fa  (Section  3.5.1). 


Figure  4,  Size  parameter  in  foam  x  (Equation  (7)  in  Table  1)  as  a  function  of  foam  void 
fraction  at  various  frequencies  and  bubble  radius  a  =  0.3  mm  at  the  peak  of  the  bubble  size 
distribution  obtained  from  oceanographic  measuremetns. 


The  figure  shows  that  x  increases  as  foam  changes  from  dry  to  wet.  This  change  of  x  in  foam  depth 
is  stronger  for  lower  frequencies  in  accord  to  the  variations  of  the  values  shown  in  Table  2.  For 
instance,  at  37  GHz  the  value  of  x  at/,  =  98%  changes  by  a  factor  of  3.8  compared  to  the  values  of  x  at 
fa  =  10%,  while  at  10.7  GHz  this  change  is  by  a  factor  of  6.6.  Table  3  summarizes  the  ranges  of 
variation  of  x  in  foam  thickness  using  the  values  of  the  predominant  bubble  size  for  each  sublayer 
(Columns  1  and  2  in  Table  3). 

Table  3.  Ranges  of  variation  of  size  parameter  (x),  and  foam  refractive  index  (Re{m}  and 
|Im{m}|)  in  foam  depth. 


Foam  Mixture 

X  (Figure  4) 

Re{wi}  (Figure  5(a)) 

llm{m}|  (Figure  5(b)) 

fa  =  98% 
a  =  1 0  mm 

0.2-10 

=1 

6  X  10 "^-3  X  10^^ 

fa  =  60% 

a  =  0.3  mm 

0.02-0.6 

3. 5-2. 5 

0.6-1 

/=io% 

a  =  0.05  mm 

0.01-0.2 

9-4.5 

2-4 

3.5.3.  Refractive  Index  of  Foam 

The  refractive  index  of  foam  m  is  obtained  with  Equation  (8)  in  Table  1.  Figure  5  presents  the  real 
and  imaginary  parts  of  m  (m'  and  m"  in  panels  a  and  b,  respectively)  as  a  function  of  frequency  at  void 
fractions  representative  for  each  of  the  three  sublayers — 98%,  60%,  and  10%.  At  a  fixed  void  fraction 
value,  the  dependence  on  frequency  is  noticeable.  But  the  variations  due  to  changes  of  the  void 
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fraction  are  stronger,  espeeially  for  the  imaginary  part.  As  expeeted,  in  dry  foam  {fa  =  98%)  m'  ~  \  and 
m"  features  the  lowest  values  of  0(10”^).  As  the  seawater  eontent  inereases  in  foam  depth,  so  does  the 
refraetion  index  (Table  3):  at^a  =  60%,  m'  =  0(3)  and  m"  =  0(1);  aty^  =  10%  m'  =  0(5  to  9)  and  m"  =  0(3). 

Figure  5.  Foam  refraetive  index  m  (Equation  (8)  in  Table  1)  as  a  funetion  of  frequeney  at 
foam  void  fractions  representing  dry  foam  {fa  =  98%),  wet  foam  {fa  =  10%)  and 
intermediate  stage  foam  {fa  =  60%)  at  fixed  seawater  temperature  (T,  =  20  °C)  and  salinity 
(5*=  34  psu):  (a)  Real  part,  Re{m};  (b)  Imaginary  part,  |Im{m}|. 
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5. 6.  Roughness  of  Foam  Layer  Interfaces 

Though  this  is  a  meehanieal,  not  dieleetrie,  property  of  the  sea  foam,  we  eonsider  the  roughness  of 
the  foam  layer  interfaees  because  it  is  useful  in  eonsidering  the  radiative  proeesses  in  foam 
(Seetion  4.2).  Strueturally,  foam  roughness  comes  about  as  eoneave  and  convex  surfaees  formed  at  the 
boundaries  from  the  bubble  eaps  protruding  into  the  air  and  indenting  into  the  water.  Roughness  due  to 
bubbles  at  the  foam  boundaries  has  not  been  presented  in  the  literature.  Here  we  make  an 
order-of-magnitude  estimate  using  two  smoothness  eriteria  (Seetion  2.3). 

We  assume  that  half  of  a  bubble  projeets  out  of  or  indents  into  a  referenee  surfaee  to  form  a  vertieal 
eharaeteristies  length  seale.  We  ean  thus  use  bubble  radii  to  represent  the  roughness  height  of  foam  at  the 
air  and  water  interfaees.  Bubble  radius  data  are  obtained  in  laboratory  and  field  experiments  [7,30,31]. 
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From  the  large-size  end  of  the  bubble  data  (a  from  1  mm  to  10  mm),  foam  rms  height  at  the  air-foam 
boundary  is  (7af=  1-3  mm.  At  the  foam-water  interfaee,  we  estimate  (7j\v  =  0.2  mm  using  the  small-size 
end  of  the  bubble  data,  radii  from  0.05  mm  to  1  mm.  For  Cre/and  A:(T(Seotion  2.3)  we  use  and  AfdXa. 
void  fraetion  of  10%  (Table  2)  for  the  air  and  water  interfaees,  respeetively. 

Figure  6  shows  all  these  estimates  together  over  the  range  of  eonsidered  frequeneies  for  6=  53°. 
Panel  (a)  represents  results  for  the  Fraunhofer  eriterion  (Seetion  2.3),  while  panel  (b)  shows  results  for 
the  k(7  eriterion.  At  the  air-foam  interfaee,  CJref  (curved  solid  line)  of  the  Fraunhofer  criterion  changes 
from  about  10  mm  at  1.4  GHz  to  0.4  mm  at  37  GHz.  Compared  to  the  rms  height  due  to  bubbles  at  this 
boundary,  (7af=  1.3  mm  (horizontal  solid  line),  the  Fraunhofer  criterion  identifies  dry  foam  as  a  smooth 
surface  for  frequencies  below  approximately  12  GHz  (vertical  solid  line).  For  higher  frequencies,  the 
dry  foam  would  be  a  rough  surface.  Approximately  the  same  frequency  limit  of  12  GHz  divides  the 
wet  foam  into  smooth  and  rough  surfaces  for  the  foam-water  interface  (dashed  lines). 

Figure  6.  Criteria  for  surface  smoothness  (Section  2.3)  at  the  air-foam  boundary  (solid 
lines)  and  foam-water  boundary  (dashed  lines):  (a)  Fraunhofer  criterion;  (b)  A:(7 criterion. 


0  5  10  15  20  25  30  35  40 

Frequency,  F  (GHz) 

According  to  the  kcj  criterion,  a  surface  is  smooth  if  k(7  <  0.2  and  very  rough  if  k(7  >  1.0  ([32]; 
p.  828).  Figure  6(b)  shows  how  the  ka  values  for  air-foam  and  foam-water  boundaries  (solid  and 
dashed  curves,  respectively)  compare  to  these  limits  (dotted  lines).  The  smooth/rough  separation  for 
the  boundaries  of  both  dry  and  wet  foam  occurs  at  around  6.5  GHz.  That  is,  according  to  the  ka 
criterion,  dry  and  wet  foam  surfaces  are  rough  for  a  wider  range  of  frequencies  than  that  established  by 
the  Fraunhofer  criterion. 


Remote  Sens.  2012,  4 


1173 


4,  Radiative  Processes  in  Foam 

Following  the  Kirchoff  s  law  for  thermal  radiation  [41] — the  emissivity  of  a  medium  at  thermal 
equilibrium  equals  its  absorptivity — ^we  know  that  in  order  to  model  the  medium  emissivity  we  need  to 
evaluate  the  absorption  losses  in  it.  We  recognize  two  possible  pathways  of  absorption  in  foam  layers 
(Section  4.5  in  [25]).  We  call  medium-intrinsic  the  absorption  of  the  radiation  propagating  along  the 
line  of  observation  due  to  the  presence  of  seawater,  a  lossy  medium.  We  call  scattering-driven  the 
absorption  caused  by  a  lengthened  path  of  propagation  of  the  radiation  along  the  line  of  observation 
due  to  the  scattering  among  densely  packed  bubbles.  The  scattering-driven  absorption  enhances  the 
medium-intrinsic  absorption. 

Note  that  we  consider  the  attenuation  caused  by  scattering-driven  absorption  different  from  the 
attenuation  caused  by  loss  of  radiation  scattered  out  of  the  line  of  observation.  Both  of  these  scattering 
effects  cause  attenuation  in  the  medium,  but  one  of  them  has  a  positive  consequence  for  our 
considerations  because  it  enhances  the  absorption  and  thus  contributes  to  the  emission  of  the  medium. 
Our  descriptive  analysis  does  not  allow  evaluation  of  the  relative  magnitude  of  these  two  scattering 
effects.  However,  heeding  the  fact  that  sea  foam  has  high  thermal  emission,  which  necessarily  derives 
from  strong  absorption  losses,  we  assume  that  the  lost  of  radiation  available  for  absorption  via 
scattering  is  small  in  foam. 

The  processes,  which  “deliver”  the  EM  radiation  available  for  absorption  in  the  foam  are  reflection 
and  surface  scattering  at  the  air-foam  and  foam-water  boundaries,  scattering  within  the  foam,  and 
transmission  through  the  foam.  In  the  following,  we  analyze  the  foam  dielectric  properties  presented  in 
Section  3  and  gain  physical  understanding  of  how  these  fundamental  radiative  processes  acquire 
unique  features  because  of  the  presence  of  vertical  void  fraction  profile  in  a  foam  layer. 

4.1.  Foam  Reflection  and  Transmission 

4.1.1.  Foam  as  Impedance  Matching 

In  absence  of  foam,  the  air  water  interface  is  highly  reflective.  The  reason  for  this  is  the  large 
difference  between  the  dielectric  constants  of  air  (eo  =  1  ^  fO)  and  seawater,  which  leads  to  very 
different  intrinsic  impedances  (Figure  2(a),  solid  and  dotted  lines).  Because  the  impedance  difference 
between  two  media  determines  the  reflectivity  at  their  interface,  r  oc  At]  (e.g.,  p.  228  in  [35]),  the 

reflectivity  at  the  air  water  boundary  is  large.  Thus,  only  a  small  portion  of  down-welling  atmospheric 
radiation  impinging  on  the  air  water  interface  would  pass  and  be  available  for  absorption  within  the 
foam-seawater  system.  As  Kirchoff  s  law  holds,  a  foam  free  water  surface  has  low  emission.  The 
cartoon  in  Figure  7(a)  illustrates  this. 

The  situation  changes  when  foam  is  present  at  the  water  air  interface  because  the  foam  acts  as 
matching  impedance.  Dry  foam  (e.g.,  fa  =  98%)  at  the  top  of  the  foam  layer  has  dielectric  properties  close 
to  those  of  air,  while  wet  foam  (fa  =  10%)  at  the  bottom  of  the  foam  layer  has  dielectric  properties  similar 
to  those  of  seawater  [25].  Accordingly,  the  intrinsic  impedances  for  wet  and  dry  foam  are  close  to  those 
of  seawater  and  air,  respectively.  Figure  2(a)  (dashed  and  dash-dotted  lines).  Therefore,  as  the  foam  void 
fraction  gradually  changes  from  the  top  of  the  foam  layer  to  its  bottom,  it  provides  a  seamless  matching 
of  the  dielectric  properties  of  air  and  seawater.  Figure  2(b).  This  action  of  a  foam  layer  on  the  water 
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surface  is  the  natural  eounterpart  of  the  teehnology  used  to  reduee  refleetion  in  optieal  systems  ([42], 
Chapter  2,  §16).  That  is,  the  foam  layer  aets  as  an  antirefleetive  eoating  on  the  sea  surfaee  allowing 
effeetive  transmission  of  down-welling  atmospherie  radiation  into  the  foam-seawater  system.  As  more  of 
the  transmitted  radiation  is  absorbed  in  the  foam-water  system,  more  is  emitted,  ultimately  resulting  in 
high  emissivity  from  foam-eovered  water  surfaees.  The  eartoon  in  Figure  7(b)  illustrates  this  situation. 


Figure  7.  Sehematie  representation  of  refleetion  and  transmission  of  ineident  (left)  and 
emitted  (right)  radiation  under  different  eonditions  at  air-sea  interfaee;  (a)  foam-free  interfaee; 
(b)  foam-eovered  interfaee;  (c)  bubbly  mixture  in  the  water  without  foam  on  the  surfaee. 


(a)  (b)  (c) 


4.1.2.  Effeetiveness  of  Impedanee  Matehing 

Williams  [19]  suggested  that  the  most  effeetive  impedanee  matehing  of  sea  foam  would  neeessarily 
involve  a  surfaee  layer.  Figure  2(b)  eonfirms  this  suggestion  by  showing  that  the  impedanee  differenee 
between  seawater  and  air  is  best  minimized  only  when  the  full  range  of  void  fraetion  is  present,  from 
~  100%  to  ~  0%.  This  eould  be  fulfdled  only  when  foam  floats  on  the  surfaee.  Strong  emission  would 
not  be  observed  when  bubble  plumes  linger  in  the  water  eolumn  but  foam  is  not  present  at  the  air-water 
interfaee  above  them  (Figure  7(e)).  The  reason  is  that  the  lower  impedanee  of  the  aerated  seawater  in 
the  bubble  plume  is  “shielded  from  the  surfaee  for  mierowaves”  ([19],  p.  223)  and  thus  its  matehing 
role  is  obviated. 

Besides  presenee  or  laek  of  foam  on  the  surfaee,  another  possibility  that  might  influenee  the 
effeetiveness  of  the  impedanee  matehing  is  that  the  vertieal  profiles  in  foam  layers  with  different 
thieknesses  may  eover  limited  instead  of  the  full  range  of  values.  For  instanee,  aetive  whiteeaps  are 
thiek  thus  they  are  antieipated  to  have  impedanee  eorresponding  to  wider  range  of  void  fraetions,  from 
an  upper  value  of  80-100%  down  to  about  10%.  Meanwhile,  the  foam  layers  of  residual  whiteeaps  are 
thinner  and  may  have  a  lower  impedanee  assoeiated  with  a  narrower/,  range  shifted  toward  lower 
values,  Q.g.,fa  e  (<60%  to  0%)  ([25],  Seetion  4.4).  Figure  2(b)  suggests,  however,  that  even  foam  with 
a  limited/,  range  eould  provide  impedanee  matehing.  It  may  not  be  the  most  effeetive,  but  still  it  eould 
mediate  some  deerease  in  the  impedanee  mismateh  between  air  and  seawater.  The  most  important 
requirement  therefore  for  the  impedanee  matehing  is  the  foam  to  float  on  the  surfaee. 
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4.2.  Surface  Scattering  of  Foam 

Only  smooth  boundary  invokes  a  specular  reflection  as  that  illustrated  in  Figure  7(a).  In  real 
geophysical  configurations,  all  boundaries  are  more  or  less  rough.  This  results  in  a  surface  scattering 
pattern,  which  adds  to  or  completely  supersedes  the  specular  reflection. 

The  strength  of  surface  scattering  is  defined  by  the  dielectric  contrast  of  two  media  on  both  sides  of 
an  interface  (as  in  Section  4.1)  and  the  degree  of  roughness  of  this  interface.  In  Section  4.1  we  showed 
that,  by  the  virtue  of  impedance  matching,  the  foam  boundaries  are  weakly  to  non-refiective  surfaces. 
This  relatively  weak  reflection  might,  however,  be  augmented  if  these  boundaries  are  rough.  In  this 
section  we  address  the  question  whether  the  roughness  of  the  foam  boundaries  is  large  enough  to  affect 
the  weak,  dielectrically  determined  surface  reflection  of  the  sea  foam. 

Figure  6  showed  that,  according  to  the  Fraunhofer  criterion,  both  dry  and  wet  foams  appear  rough  to 
EM  radiation  above  12  GHz,  while  the  kcr  criterion  places  this  limit  to  about  6.5  GHz.  To  put  into 
prospective  these  results,  we  compare  the  foam  rms  height  values  <7af  and  Ofiv  used  in  Figure  6  to  the 
rms  height  of  a  wind-roughened  sea  surface  (7u.  We  obtain  CTu  from  the  variance  of  the  sea  surface 
displacement  h  about  the  mean  surface  at  a  given  point.  It  is  defined  as  the  zeroth  moment  of  the  wave 

spectrum  S{K),  (T^=m^=  | S(K)dK  =  {h^'^  [43,44],  where  K  =  2nlA  is  the  wave  number  of  an  ocean 

^low 

wave  with  wavelength  A.  Note  that  CJu  is  a  measure  of  the  geometric  roughness  of  a  wind-driven  sea 
surface,  not  the  aerodynamic  roughness  usually  represented  by  the  roughness  length  zo  [45,46]. 

The  choice  of  integration  limits  Kiow  and  Kup  above  determines  the  wavelength  scales  contributing  to 
(7u.  But  these  length  scales,  determined  mostly  by  geometric  consideration,  also  need  to  be  tied  to  length 
scales  important  for  the  aerodynamic  surface  roughness  because  this  is  the  variable  used  for 
air-sea  interaction  studies.  Our  reasoning  in  choosing  Kiow  and  Kup  is  thus  as  follows.  First,  note  that 
remote  sensing  of  the  sea  surface  is  effective  when  the  probing  EM  wavelengths  Ao  are  comparable  to  the 
length  scales  of  the  ocean  waves.  Eor  example,  for  scatterometers,  ocean  length  scales  ranging  from 
about  40  cm  to  less  than  2  cm  correspond  to  Ao  from  25  cm  to  3  cm  at  incidence  angles  of  20°  and  60° 
([36],  p.  1705).  On  this  premiss,  one  can  expect  the  probing  wavelengths  Ao  for  the  range  of  frequencies 
investigated  here  (Table  2)  to  be  most  relevant  to  ocean  lengths  scales  from  about  40  cm  to  a  few 
millimeters.  These  length  scales  represent  short  gravity  and  capillary  water  waves.  Next,  note  that 
analysis  of  the  mean  square  slope  of  wave  number  spectrum,  resulting  from  oceanographic  surface  wave 
measurements,  has  shown  that  intermediate-scale  waves  with  A  between  6  m  and  2  cm  are  the  dominant 
contributors  to  the  aerodynamic  surface  roughness  in  oceans  [47].  Combining  these  two  considerations, 
we  restrict  our  estimates  to  ocean  length  scales  from  40  cm  to  2  cm.  The  corresponding  wavenumbers  K 
range  from  3.14  rad-cm~^  to  0.16  rad-cm~\  To  cover  this  range  of  K  values,  we  obtained  mo  for  three 
portions  of  the  wave  spectrum;  S2{IC),  S^iK),  and  SfK)  as  defined  in  [43].  In  the  estimates  we  used  a 
friction  velocity  u*  of  24  cm-s~'  (Uio  ~  7  m-s~^),  a  value  close  to  the  globally  averaged  wind. 

Eor  the  considered  range  of  length  scales  for  ocean  waves,  we  find  that  the  wind-roughened  rms 
height  ffu  ranges  from  about  1  mm  for  the  short  end  of  water  waves  (A  of  2  cm)  to  5.4  mm  for  longer 
water  waves  (A  of  40  cm).  We  obtain  similar  estimates  for  if  an  empirically  formulated  wave 
spectrum  in  terms  of  mean  square  slope  is  used  [47]  instead  of  S{IC). 
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A  comparison  of  these  estimates  to  <7a/=  1.3  mm  implies  that  there  eould  be  eases  in  whieh  the 
roughness  of  the  dry  foam  eould  be  eomparable  to  the  roughness  ereated  by  the  shortest  eapillary 
waves.  But  the  dieleetrie  eontrast  between  dry  foam  and  air  is  almost  none  existing  (Figure  2  and 
Seetion  4.1).  This  leads  ultimately  to  laek  of  surfaee  seattering  at  the  air-foam  boundary. 

For  the  wet  foam  boundary,  its  roughness  {Ofiv  =  0.2  mm)  is  an  order  of  magnitude  below  that  of  the 
shortest  eapillary  waves  (1  mm).  This  implies  that  wet  foam,  whether  at  the  surfaee  or  as  a  sublayer  in 
the  foam  thiekness,  is  usually  smoother  than  the  sea  surfaee  roughness.  We  thus  deem  the  wet  foam 
boundary  as  ineapable  of  providing  signifieant  surfaee  seattering  as  well. 

Our  eonelusion  of  negligible  surfaee  seattering  by  any  of  the  foam  layer  boundaries  is  eonsistent 
with  the  general  observation  that  the  refleetivity  of  a  foam-eovered  surfaee  is  not  eomparable  to  that  of 
foam- free  rough  sea.  In  faet,  mierowave  models  (e.g.,  [48])  usually  assume  that  “foam  tends  to 
deerease  the  refleetivity”  ([49],  p.  24).  Note  that  the  eonelusions  in  this  seetion  result  from  the 
eombined  effeet  of  the  geometrie  roughness  and  the  dieleetrie  eontrasts  at  the  boundaries  of  a  foam 
layer  with  well-stratified  void  fraetion.  Deviations  from  these  eonelusions  are  possible  when  the  void 
fraetion  profile  is  limited  (Seetion  5.2). 

4.3.  Weak  Volume  Scattering  Throughout  a  Foam  Layer 

While  we  reeognize  the  presenee  of  seattering-driven  pathway  of  attenuation  in  foam  and  expeet  it 
to  play  an  important  role  in  a  strueture  eondueive  to  signifieant  seattering  (namely,  the  elose  paeking 
of  bubbles),  a  previous  survey  ([24],  Seetion  2.3)  has  shown  that  at  the  eonsidered  frequeneies  the 
seattering  within  the  sea  foam  is  weak.  Here,  we  explain  this  eounterintuitive  eonelusion. 

Generally,  one  eonjeetures  on  the  presenee,  or  laek  of,  volume  seattering  in  a  medium  by  eonsidering 
its  effeetive  penetration  depth  and  the  spatial  and  dieleetrie  irregularities  in  its  volume  ([32],  Chapter  1 1^). 
The  spatial  and  dieleetrie  irregularities  are  usually  quantified  with  their  spatial  distribution  and  dieleetrie 
eontrast.  The  size  parameter  and  the  refraetive  index  of  a  medium  express  these  two  quantities.  Thus,  in 
the  following  we  analyze  S,  x  and  m  for  elues  regarding  the  volume  seattering  in  foam. 

Our  reasoning  proeeeds  in  two  steps.  First,  employing  general  seattering  rules,  we  establish  the 
expeeted  seattering  behavior  for  various  foam  struetures.  Seeond,  we  show  how  this  expeeted  behavior 
is  altered  (or  not  followed)  in  natural,  vertieally  stratified  foam  layers. 

4.3.1.  Expeeted  Seattering  in  Foam  Layers 

Using  S,  we  ean  identify  the  strueture  that  eould  eontribute  to  volume  seattering  the  most.  Beeause 
eompared  to  fixed values  the  integral  effeet  of  a  vertieaiy^  profile  is  to  limit  d  [25],  we  infer  from 
Figure  1  that  signifieant  volume  seattering  eould  be  expeeted  in  aetive  whiteeaps  as  they  involve  more 
dry  foam  and  thus  have  the  largest  penetration  depth.  Beeause  the  penetration  depth  of  intermediate-wet 
and  wet  foam  {fa  <  60%,  e.g.,  residual  whiteeaps)  is  mueh  smaller,  the  inhomogeneities  within  eould  be 
assumed  to  eontribute  little  to  volume  seattering. 

Using  X  and  m,  we  ean  further  assess  the  strength  and  relative  importanee  of  the  volume  seattering  in 
dry  and  wet  foam.  For  this  purpose,  we  employ  two  basie  rules  established  for  seattering  by  partieles. 
Namely,  (i)  the  larger  the  partiele  eompared  to  the  wavelength  {^«  a  ox  x»  1)  the  more  effeetive  the 
seattering  ([37],  p.  68);  and  (ii)  the  larger  the  dieleetrie  eontrast  between  the  refraetive  indiees  of  the 
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particle  and  the  surrounding  medium  the  stronger  the  seattering.  That  is,  in  the  ease  of  partieles  in  air, 
there  is  no  scattering  for  m  =  1  and  scattering  is  small  for  m  near  1  ([37],  p.  172).  We  apply  these  rules 
for  partieles  to  foam  on  the  basis  of  the  connection  between  foam  layers  and  bubbles  (Seetion  3.5.1). 

We  start  with  the  analysis  of  x  values.  Figures  8a-c  show  x  as  a  function  of  frequency  in  dry, 
intermediate-wet,  and  wet  foams  for  bubble  sizes  characteristic  for  each  sublayer  (Figure  3  and  Table  3). 
In  eaeh  panel,  the  size  parameter  assoeiated  with  the  eharaeteristie  bubble  size,  whieh  is  the  least  likely 
for  the  depicted  void  fraction  is  plotted  in  gray.  We  observe  two  trends  from  panels  (a)  to  (c)  in 
Figure  8.  One  trend  is  that  for  most  frequeneies,  foam  structures  eomp rising  large  bubbles  (a  =  10  mm) 
have  X  >  1  or  X  »  1  (dashed  lines),  whatever  the  void  fraction  those  bubbles  form  by  having  various 
wall  thicknesses  w  (reeall  Section  3.5.1).  According  to  rule  (i)  above,  this  means  that  foam  containing 
larger  bubbles  would  have  strong  scattering.  Another  trend  observed  in  Figure  8  is  that  x  »  1  in  wetter 
foam  (panels  (b)  and  (c))  than  in  dry  foam  (panel  (a)).  From  these  two  trends,  we  surmise  that  significant 
and  effective  volume  seattering  would  be  present  in  wet  foam  comprising  large,  thiek-walled  bubbles; 
Figure  8(c)  (dashed  line)  exemplifies  this  generalization. 

4.3.2.  Altered  Scattering  in  Foam  Layers 

The  expected  scattering  behavior  established  in  Section  4.3.1  is  not  readily  applicable  for  natural 
foams  in  which  bubbles  sort  themselves  by  size  in  the  foam  depth  (Section  2.1).  Because  of  the  bubble 
size  stratifieation,  large  bubbles  are  the  least  likely  to  oecur  or  non-existing  in  wetter  sublayers; 
Figure  8(c)  shows  this  (the  dashed  line  is  gray).  Thus,  effective  volume  scattering  should  be  confined 
to  mostly  dry  foam.  However,  dry  foam  has  the  lowest  x  values  (Figure  8(a))  eompared  to  x  of  other 
foam  structures  (Figure  8(b,c)).  Analysis  of  x  thus  reveals  the  first  limitation  for  achieving  strong 
volume  seattering  in  natural  foam  layers.  Namely,  the  volume  seattering  is  restricted  to  the  least 
effective  seattering  structure,  the  dry  foam. 

Analysis  of  m  reveals  another  limitation  for  effective  volume  scattering  in  natural  foam  layers. 
From  Figure  5  and  rule  (ii)  (Seetion  4.3.1),  we  see  that  the  x-limited  effectiveness  of  the  dry  foam  as  a 
scattering  medium  is  further  diminished  by  its  small  dielectric  contrast.  Having  m  ~  1  (solid  line  in 
Figure  5),  the  refraetive  index  of  the  dry  foam  does  not  differ  mueh  from  that  of  the  air.  Thus,  the 
scattering  in  dry  foam  is  negligible.  Meanwhile,  as  m  inereases  in  foam  depth,  the  perceived 
ineffeetiveness  of  the  wet  foam  as  a  scattering  medium  based  on  x  «  1  could  be  somewhat 
strengthened.  This  is  espeeially  true  when  wet  foam  is  adjaeent  to  air  (residual  whitecaps  might 
exemplify  such  a  case).  But  when  wet  foam  is  part  of  a  vertically  structured  foam  layer  (as  in  active 
whiteeaps),  the  inereased  values  of  m  would  still  lead  to  modest  enhaneement  of  the  volume  seattering 
because  the  immediately  adjacent  sublayers  have  close  m  values. 

Overall,  we  infer  that,  in  a  foam  mixture  eomprising  a  distribution  of  bubble  sizes  sorted  in  foam 
depth,  the  simultaneous  decrease  of  and  the  number  of  bubbles  with  radii  effeetive  for  seattering  keeps 
the  X  values  relatively  low  and  confined  in  a  narrow  range.  As  a  result,  the  scattering  in  foam  remains 
less  effeetive.  This  is  eompounded  with  the  effect  of  the  dielectric  contrast,  which  is  expressed  via  m 
increasing  in  foam  depth  yet  staying  adjacent  to  sublayers  with  similar  m  values.  The  scattering  in  wet 
foam  is  somewhat  more  effective  than  the  scattering  in  dry  foam,  however  the  thiekness,  thus  the 
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penetration  depth,  of  wet  foam  restricts  this  further.  The  overall  outcome  is  weak  scattering  throughout 
the  foam. 


Figure  8,  Size  parameter  in  foam  x  (Equation  (7)  in  Table  1)  as  a  function  of  frequency 
and  foam  void  fraction  representing;  (a)  dry  foam  {fa  =  98%);  (b)  intermediate-wet  foam 
{fa  =  60%)  and  (c)  wet  foam  {fa  =  10%).  For  each  /«  region  the  size  parameters  for  three 
possible  bubble  radii  are  shown.  The  radii  least  likely  for  each  type  of  foam  are  in  gray. 
Fixed  seawater  temperature  (T,  =  20  °C)  and  salinity  {S  =  34  psu). 
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4.4.  Strong  Absorption  by  Wet  Foam  Revisited 

Section  4.3  provides  a  physical  explanation  of  why,  of  the  two  pathways  for  absorption  losses  in 
foam,  the  scattering-driven  one,  though  present,  is  not  a  major  contributor.  This  thus  qualitatively 
confirms  the  conclusion  Anguelova  and  Gaiser  ([25];  Section  4.5)  reached  that  the  medium-intrinsic 
absorption  is  the  defining  controlling  factor  for  high  foam  absoprtivity.  Anguelova  and  Gaiser  [25] 
also  showed  that  the  thin  water-laden  wet  foam  could  provide  a  strong  emission  from  foam-covered 
sea  surfaces.  Here  we  give  further  evidence  for  this  conclusion  using  the  results  for  x  and  m. 

Besides  small  skin  depth  ([25];  Section  4.1),  another  quantitative  measure  for  the  high  absorptivity 
of  foam  as  a  medium  is  the  fact  that  the  refractive  index  (and  permittivity  ^  )  of  foam  is  a  complex 
number  (m  =  m' -im" ,  Section  3.5.3).  Media  with  even  very  small  m"  values  have  been  called  strongly 
absorbing.  On  the  basis  of  m"  ^  0  and  m  -  1  «  1,  which  means  that  a  negligible  fraction  of  the 
radiation  is  reflected,  one  could  conjecture  that  even  dry  foam  (Figure  5,  solid  lines)  would  be 
absorptive.  In  fact,  a  combination  of  refractive  index  values  similar  to  that  of  dry  foam  with  large  size 
parameter  (e.g.,  x  >  100)  could  represent  a  medium  whose  behavior  approaches  that  of  a  black  body 
([37],  p.  269). 

However,  such  a  combination  of  m  and  x  is  not  observed  in  natural  dry  foam  because  x  does  not 
exceed  10  (Figure  8(a)).  The  combination  of  a  narrow  range  of  x  (x  <  10)  with  m"  «  1  thus  excludes 
the  possibility  that  dry  foam  contributes  to  medium-intrinsic  absorption  losses.  Meanwhile,  m'  and  m" 
values  for  wet  foam  are  of  the  same  order  (Figure  5,  dashed  lines)  signifying  that  this  is  an  intrinsically 
absorptive  medium.  However  little  the  volume  scattering  in  wet  foam  (Section  4.3),  it  adds  to  the 
medium-intrinsic  absorption,  making  the  wet  foam  a  more  effective  absorber  than  the  dry  foam. 

5,  Discussion 

5.1.  Concept  for  the  High  Foam  Emissivity 

Various  suggestions  can  be  found  in  the  literature  explaining  the  high  emissivity  of  foam-covered 
surfaces.  Williams  [50]  was  the  first  to  suggest  the  importance  of  the  impedance  matching  that  the  sea 
foam  provides  between  the  air  and  seawater.  Consequent  experiments  led  to  the  tentative  conclusion 
that  roughness  expressed  as  “distortion  of  the  water  surface  at  the  interstices  between  the  bubbles  at 
the  bottom  of  the  foam  layer”  ([51],  p.  3097)  together  with  the  dielectric  constant  of  foam  “produce  a 
zone  of  radically  tapered  dielectric  constant”  ([19],  p.  223;  cited  also  by  [21];  [36],  p.  1454;  [48]).  The 
Russian  literature  on  the  subject  approaches  the  problem  from  a  microscopic  point  of  view  and 
promotes  the  explanation  that  the  bubbles  within  the  foam  “represent  some  kind  of  small  black  bodies, 
which  intensively  absorb  the  EM  energy  owing  to  inter-bubble  diffraction”  ([52],  p.  511).  Describing 
the  sea  foam  as  “a  nonreflecting  ...  absorbing  but  partially  transparent  medium”,  Wilheit  ([48],  p.  246) 
touched  upon  the  effective  transmission  of  EM  radiation  through  foam.  Williams  [19]  measurements 
involving  foam  on  different  underlying  substrates,  water  or  metal,  suggest  that  not  just  the  foam  itself, 
but  also  the  system  of  foam  floating  on  seawater  is  the  one  that  could  lead  to  high  emissivity. 

This  brief  account  of  proposed  explanations  demonstrates  that  one  or  another  element  of  the  physical 
system  that  a  foam-covered  sea  surface  represents  has  been  identified  and  used  to  interpret  the  observed 
high  emissivity.  None  of  these  explanations,  however,  combines  these  elements  into  a  consistent  whole. 
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Meanwhile,  the  analysis  in  this  study  of  the  radiative  proeesses  in  foam  (Seetions  4. 1-4.4)  shows  that 
none  of  these  separate  elements  alone  is  enough  to  aeeount  for  the  high  emissivity  e/of  sea  foam.  Rather, 
all  foam  properties  and  all  radiative  proeesses  engage  in  an  emergent  behavior  {i.e.,  the  whole  is  greater 
than  the  sum  of  its  parts  [53])  and  eontribute  in  some  way  to  bring  about  the  overall  effeet  of  a  high  e/. 
Moreover,  we  argue  that  what  makes  the  sea  foam  an  exeeptional  emitter  is  the  eombination  of  these 
properties  and  proeesses  in  a  vertieally  struetured  foam  layer  floating  on  seawater  surfaees.  The 
following  is  a  summary  of  our  eoneeptual  understanding  of  how  this  happens. 

First,  two  eomponents  ensure  the  high  emissivity  of  foam-eovered  oeean  surfaees:  (1)  High 
absorption  losses  in  a  eoupled  foam-seawater  system;  and  (2)  effeetive  transmission  of  EM  radiation 
through  the  foam.  If  one  of  these  eomponents  were  missing,  a  high  ef  eould  not  be  guaranteed.  The 
potential  of  the  foam-seawater  system  for  strong  absorption  and  emission  eannot  be  realized  without 
effeetive  transmission  and  propagation  of  radiation  aeross  layer  boundaries  into  the  foam;  and  vice 
versa — effeetive  transfer  aeross  the  layer  boundaries  is  not  enough  if  a  highly  absorptive  agent  is  not 
in  plaee  to  earry  on  high  absorption  and  emission. 

Seeond,  absorption  losses  in  foam  are  provided  by  medium-intrinsie  and  seattering-driven  absorption. 
Absorption  in  wet  foam  is  relatively  more  important  than  absorption  in  dry  foam.  Furthermore,  though 
strong,  the  absorption  of  the  foam  layer  itself  might  remain  somewhat  limited.  The  seawater  below  the 
foam  layer  is  a  major  absorptive  agent  whose  emission  adds  to  that  of  the  foam  layer  to  yield  the 
highest  absorption  losses  in  a  eoupled  foam-seawater  system. 

Third,  with  its  intrinsie  impedanee,  foam  on  the  oeean  surfaee,  in  any  stage  of  its  lifetime,  provides 
impedanee  matehing  between  the  eontrasting  dieleetrie  properties  of  air  and  seawater.  This  deereases 
refleetions  at  the  air-foam  and  foam-seawater  boundaries  and  assists  effeetive  entry  of  EM  radiation 
into  the  foam-seawater  system.  Existing,  yet  weak  surfaee  seattering  at  these  boundaries  and  weak 
volume  seattering  within  the  foam  provide  the  maximal  portion  of  the  propagating  radiation  to  be 
delivered  to  and  absorbed  by  the  wet  foam  and  seawater. 

Forth,  the  foam-seawater  system  reaehes  its  maximal  effeetiveness  and  beeomes  a  blaek-body-like 
emitter  when  a  vertieal  profile  of  properties  in  the  depth  of  a  foam  layer  eontrols  the  relative  importanee 
and  eontribution  of  the  radiative  proeesses.  Only  gradual  inerease  of  seawater  eontent  and  gradual 
ehange  of  spatial  inhomogeneities  minimize  radiation  attenuation  in  dry  foam  and  maximize  the 
absorption  losses  in  wet  foam  and  underlying  seawater.  Only  a  gradual  ehange  of  the  foam  intrinsie 
impedanee  over  the  widest  possible  range  of  values,  from  that  of  the  air  to  that  of  the  seawater,  provides 
weakly  refieetive  to  non-refleetive  foam  boundaries  and  maximizes  the  entry  of  EM  radiation. 

In  short,  only  in  vertieally  struetured  foam  layers  every  property  and  proeess  eonspires  to  bring 
maximum  EM  radiation  to  the  most  effeetive  absorptive  agent  and  thus  produeing  maximum  absorption 
and  emission. 

5.2.  Foam  as  a  Dynamic  System 

Our  diseussion  of  the  sub-surfaee  bubble  plume  and  the  overlying  foam  layer  generated  is  presented 
up  to  this  point  as  a  statie,  time-independent  system.  The  impaet  that  the  dynamie  of  the  system  may 
have  on  the  foam  strueture  and  its  emissivity  was  alluded  to  in  our  eonsiderations  of  the  effeetiveness 
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of  the  impedance  matching  provided  by  whitecaps  in  their  active  and  residual  lifetime  stages 
(Section  4.1.2).  Here  we  address  further  the  temporal  variations  of  the  system. 

The  rational  for  these  considerations  is  that  the  void  fraction  in  the  bubble  clouds  formed  after  wave 
breaking  varies  by  roughly  two  orders  of  magnitude  between  initial  entrainment  and  the  loss  of  large 
bubbles,  an  interval  of  roughly  a  second  or  so.  Indeed,  in  a  laboratory  study,  Anguelova  and  Huq  [59] 
observed  disintegration  of  the  initial  compact  form  of  a  whitecap  into  winding  streaks  of  foam  after 
only  100  ms.  Such  a  time  scale  is  comparable  to  the  time  period  of  90  ms  reported  by  Deane  and 
Stokes  [60]  for  the  dispersion  of  a  compact  mass  of  bubbles  into  individual  ones  in  a  surf  zone.  For 
about  one  wave  period  (about  1  s),  Anguelova  and  Huq  [59]  documented  a  decrease  of  the  void 
fraction  from  near  unity  to  about  20%.  They  also  observed  that  the  evolution  of  the  bubble  plume,  and 
the  foam  layer  above,  could  persist  for  many  wave  periods. 

The  question  is;  What  is  the  impact  of  such  dynamical  changes  on  a  field  of  whitecaps  in  the  open 
ocean?  Statistical  analysis  of  void  fraction  values  representative  for  the  initial  air  entrainment  and  for  the 
subsequent  decaying  stage  suggest  the  presence  of  both  thick  foam  layers  with  high  void  fractions  of 
70%-90%,  and  thinner,  wetter  foam  layer  with  lower  void  fraction  values  (<80%)  (e.g..  Figures  12  and  13 
in  Anguelova  and  Huq  [59]).  Because  evolving  bubble  clouds  keep  replenishing  the  foam  layers  on  the 
surface,  their  statistics  reveal  that  the  temporal  variations  lead  to  a  high  variability  of  the  foam  layer 
properties.  That  is,  in  the  open  ocean  at  any  moment  foam  layers  with  a  distribution  of  void  fraction 
values  at  their  boundaries  could  exist.  These  would  include  combinations  of  dry  foam  next  to  a 
foam- water  boundary  or  wet  foam  next  to  air.  Therefore,  some  of  the  various  foam  structures  may  exhibit 
a  distinct  impedance  discontinuity  at  the  mean  ocean  surface;  others  may  represent  stratified  foam  well. 
This  in  turn  has  important  implications  for  the  reflection  and  transmission  of  microwave  energy. 

The  next  logical  question  is:  How  can  we  model  such  a  variability  of  the  foam  layer  properties? 
One  simple  approach  is  to  assume  that  the  foam  structure  analyzed  in  the  present  study  represents  a 
likely  scenario  in  the  early  stage  of  foam  development.  A  significant  improvement  upon  such  a  simple 
approach  would  be  to  develop  an  elaborate  physical  model  for  the  time -varying  properties  of 
sub-surface  bubble  clouds,  about  which  relatively  little  is  known,  or  the  surface  foam  they  generate, 
about  which  even  less  is  known.  A  starting  point  could  be  to  generate  a  void  fraction  profile  likely  to 
be  representative  of  the  bubble  plume/foam  system  during  whitecap  decay.  Anguelova  and 
Gaiser  ([25],  Section  4.4)  discussed  the  modeling  of  void  fraction  profiles  by  noting  that  one  needs  to 
consider  the  profile  shape,  profile  range,  and  the  profile  variations  with  foam  thickness,  and  how  these 
choices  affect  the  foam  emissivity. 

5.3.  Modeling  Scattering  in  Foam 

Foam  emissivity  can  be  modeled  following  two  approaches  [36,52,54].  The  intensity  approach 
employs  the  radiative  transfer  theory  to  calculate  the  microwave  radiance  of  a  medium.  The  wave 
approach  employs  solutions  of  Maxwell  equations  to  compute  the  scattering  losses  in  a  medium  and  use 
those  to  obtain  emissivity.  The  PS  formula  represents  the  most  rigorous  starting  point  when  obtaining  the 
scattering  of  sea  foam  is  pursued.  The  reason  is  that  the  PS  formula  is  identical  to  the  low  frequency  limit 
form  of  the  effective  permittivity  obtained  with  the  strong  permittivity  fluctuations  theory  for  scattering 
from  a  random  medium  ([54],  pp.  378-379).  According  to  Tsang  et  al.  ([54],  p.  429),  “the  applicability 
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of  the  mixing  formula  to  a  particular  dense  medium  in  the  very  low  frequency  limit  will  help  in  deciding 
the  multiple  scattering  theory  to  be  employed  at  higher  frequencies.”  Because  we  use  the  PS  mixing  rule 
in  this  study  (Section  3.1),  we  are  in  good  position  to  qualitatively  discuss  which  scattering  theory  would 
be  appropriate  when  the  wave  approach  is  chosen  to  obtain  sea  foam  emissivity. 

5.3.1.  Scattering  Regimes  in  Foam 

Aden  and  Keker  [55]  and  Dombrovskiy  [56]  extended  the  classical  Mie  theory  to  apply  it  to  hollow 
spherical  particles  (bubbles  in  our  case).  Dombrovskiy  [57]  used  this  theory  to  calculate  the  absorption 
and  scattering  of  foam  at  several  microwave  frequencies.  He  made  the  same  calculations  also  using  the 
Rayleigh  approximation  and  found  that  the  results  differ  from  those  of  the  Mie  theory  by  3%  for 
thick-walled  bubbles  at  37  GHz  and  by  no  more  than  0.3%  for  thin-walled  bubbles  at  37  GHz  and 
10  GHz.  Analyzing  the  scattering  parameters  of  foam  mixtures  (Section  3.5  and  Table  3)  instead  of 
individual  bubbles  [57],  we  also  find  that  the  Rayleight  approximation  represents  well  the  scattering 
regime  in  vertically  structured  foam  layers.  Here  we  show  how  we  reach  this  result.  In  our 
considerations,  we  use  the  so-called  m-x  domain,  which  summarizes  all  possible  scattering  regimes  ([37], 
p.  132)  depending  on  the  values  of  the  medium  size  parameter,  refractive  index,  and  their  product 
(Appendix). 

Following  conditions  (A.1-A.3),  we  evaluate  the  entries  in  Table  3  as  small  (5),  large  (/),  or  arbitrary  {arb) 
values  and  formulate  the  entries  in  Table  4  forx,  m  -  1,  andx(m  -  1).  We  compare  these  to  the  regimes  listed 
in  Table  A  and  infer  the  following  scattering  regimes  for  the  most  realistic  situations  encountered  in  foam. 

Table  4,  Scattering  regimes  for  foam  evaluated  by  comparison  with  van  de  Hulst  ([37], 

Table  8).  Vales  marked  as  s  (small),  /  (large),  or  arb  (arbitrary)  following  the  van  de  Hulst 
inequalities  ([37],  p.  132). 


Foam  sublayer 

X 

m  -  1 
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Regime 
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a  =  0.3  mm 

s 

arb 

s 

R 

Rayleigh 
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a  =  0.05  mm 

s 

arb 

s 

R 

Rayleigh 

First,  considering  that  according  to  Table  4  x{m  -  1)  is  small,  the  possible  scattering  regimes  in 
foam  are  those  in  the  upper  left  side  of  the  m  -  x  square  (left  and  above  the  dashed  line  in  Figure  9) 
excluding  the  regimes  of  total  reflection  and  geometric  optics. 

Next,  according  to  Tables  3  and  4,  x  is  small  (usually  «  1)  for  intermediate -wet  and  wet  foam,  and 
when  it  is  >  1  (in  dry  foam)  it  is  no  higher  than  10;  this  is  way  below  the  possible  upper  limit  of  x  >  50 
at  microwave  frequencies  (Appendix).  This  keeps  the  scattering  in  foam  mostly  in  the  left  part  of  the 
m-x  square.  With  very  large  x  excluded,  we  rule  out  the  anomalous  diffraction  ([37],  Section  1 1.3)  as 
a  possible  scattering  regime  in  foam. 

Finally,  values  m  -  1  in  Table  4  vary  in  foam  from  small  to  arbitrary.  Large  m  values  (up  to  9,  the 
possible  upper  limit  at  microwave  frequencies  as  noted  in  the  Appendix)  combined  with  small  x  for  the 
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intermediate-wet  and  wet  foams  lead  to  optieal  resonanee  {i.e.,  Mie  seattering),  Figure  9.  Indeed, 
strong  resonanee  peaks  in  the  extinetion  is  doeumented  for  m  =  8.9  -  /0.69  and  x  <  0.5  ([37],  Figure  60). 
But  even  a  small  inerease  of  |Im{m}|  (e.g.,  from  a  value  of  0.69  to  1.96)  effeetively  damps  any 
resonanee  effeets.  The  large  |Im{m}|  values  in  intermediate -wet  and  wet  foam  (Table  3)  thus  restriet 
the  oeeurrenee  of  elassieal  Mie  seattering  in  foam. 

Overall,  we  surmise  that  with  both  x  and  x{m  -  1)  small  and  m  -  \  restrieted  to  less  than  8,  the 
prevailing  seattering  regimes  in  foam  for  the  range  of  eonsidered  frequeneies  are  Rayleigh  and 
Rayleigh-Gans  (also  ealled  Rayleigh-Debye  or  Born  approximation)  regimes.  In  the  m  -  x  domain 
(Figure  9),  this  is  roughly  the  region  to  the  left  of  the  dotted  line  and  above  the  dash-dotted  line. 


Figure  9.  Seattering  regimes  (aeronyms  in  Table  A)  in  the  m  -  x  domain.  A  simplified 
version  of  the  m  -  x  domain  given  in  van  de  Hulst  ([37],  Figure  20). 
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5.3.2.  Variations  of  the  Seattering  Regimes  in  Foam 


Anything  that  ehanges  the  values  of  the  seattering  parameters  (x,  m,  and  x{m  -  1))  in  foam,  would 
eause  deviations  from  the  eonelusion  in  Seetion  5.3.1.  Of  eourse,  the  foremost  reason  of  seattering 
regime  ehange  in  a  foam  layer  is  its  vertieal  strueture.  Going  from  the  foam-air  interfaee  toward  the 
foam-water  boundary,  the  seattering  regime  in  foam  would  ehange  from  Rayleigh-Gans  seattering  in 
dry  foam  (80%  <  fa  ^  100%)  to  Rayleigh  seattering  in  wet  foam  (0%  <  fa  ^  40%)  with  a  gradual 
transition  between  these  two  in  the  intermediate-wet  foam  (40%  ^fa'^  80%).  Other  possible  reasons 
that  eould  ehange  the  values  of  the  seattering  parameters  are;  (i)  Changes  in  m  due  to  full  aeeount  of 
the  seattering  in  e/;  (ii)  ehanges  inx  due  to  SST  effeet  (at  fixed  salinity)  on  ^and  with  this  on  Af;  and 
(iii)  ehanges  in  x  due  to  salinity  effeet  (at  fixed  SST)  on  a  and  ^(thus  Af. 

To  evaluate  the  expeeted  ehanges  in  the  seattering  regimes  due  to  (i),  let  us  assume  that  by  ignoring 
seattering  we  deviate  from  the  total  attenuation  (absorption  +  seattering)  ke  by  overestimating  or 
underestimating  it.  Beeause  of  the  direet  proportionality  Im(.^^^)  =  m"  (Equations  (2),  (3),  and  (8) 


in  Table  1),  m"  of  intermediate-wet  and  wet  foams  (Table  3)  would  inerease  or  deerease.  Meanwhile 
the  values  of  Re{m}  are  expeeted  to  ehange  little  eompared  to  those  reported  in  Table  3. 

All  values  used  in  the  eonsiderations  until  now  are  for  SST  of  20  °C  and  salinity  of  34  psu.  SST 
does  not  influenee  the  bubble  size  at  fixed  salinity.  However,  the  temperature  dependenee  of  Sf  may 
inerease  notieeably  at  lower  SST  ([24],  Figure  6).  Sueh  a  result  would  deerease  size  parameter  x, 
whieh  would  eonfine  the  seattering  in  foam  to  the  transitional  (R-RG  in  Figure  9)  or  Rayleigh  regimes. 
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Salinity  affects  both  Xf  (via  the  salinity  dependence  of  ej)  and  the  bubble  size  a.  Seawater  with 
lower  salinity  (at  fixed  SST)  breaks  into  fewer  but  larger  bubbles  (e.g.,  [58],  Figure  5).  Meanwhile, 
Anguelova  ([24],  Figure  8)  shows  that,  as  salinity  decreases,  the  real  part  of  foam  permittivity 
increases,  negligibly  so  for  higher  frequencies  (18  GHz  to  37  GHz)  and  more  noticeably  for  lower 
(1.4  GHz  to  10  GHz).  Both  changes — shift  of  bubble  size  distribution  toward  larger  sizes  and  shorter 
wavelengths — lead  to  an  increasing  x  at  lower  salinity.  In  addition,  lower  salinity  values  decrease 
seawater  conductivity  thus  the  imaginary  part  of  foam  permittivity  and  refraction  index  m.  The  real 
parts  of  ^and  m  would,  however,  keep  the  m  -  1  values  no  higher  than  about  10.  Therefore,  at  lower 
salinity  the  scattering  regime  would  shift  more  toward  the  Rayleigh-Gans  (Born  approximation) 
regime.  Similar  reasoning  points  out  that  higher  salinity  reinforces  conditions  for  predominantly 
Raylegh  scattering  in  foam. 

6.  Conclusions 

Passive  microwave  remote  sensing  of  whitecap  fraction  at  frequencies  from  1  to  37  GHz  requires 
modeling  of  ocean  surface  emissivity  fully  covered  with  foam  e/.  To  model  e/,  knowledge  of 
mechanical  and  dielectric  properties  of  sea  foam  is  necessary.  Here  we  present  a  full  physical 
description  of  the  dielectric  properties  of  sea  foam  layers  with  a  vertical  profile  of  void  fraction.  These 
include  foam  dielectric  constant,  foam  skin  depth,  foam  impedance,  wavelength  variations  in  foam 
depth,  roughness  of  foam  layer  interfaces  with  air  and  seawater,  and  foam  scattering  parameters  such 
as  size  parameter,  and  refraction  index. 

Foam  dielectric  properties  are  analyzed  to  gain  physical  insights  about  scattering,  absorption, 
reflection  and  transmission  in  foam  (Section  4).  These  fundamental  radiative  processes  acquire  unique 
features  in  vertically  stratified  foam  layers,  namely: 

(1)  Weak  or  negligible  reflection  at  the  foam  layer  interfaces  with  air  and  seawater  due  to  foam 
impedance  matching. 

(2)  Floating  foam  layers  provide  the  most  effective  impedance  matching. 

(3)  Weak  or  negligible  surface  scattering  at  the  foam  layer  interfaces  due  to  small  roughness  and 
low  dielectric  contrast  of  foam  boundaries. 

(4)  Weak  volume  scattering  throughout  the  foam  due  to  simultaneous  decrease  in  the  foam  depth  of 
radiation  wavelength  and  bubbles’  radii  effective  for  scattering. 

(5)  Absorption  losses  are  predominantly  confined  to  the  wet  portion  of  the  foam  layer. 

We  put  all  these  elements  together  and  offer  a  conceptual  understanding  for  the  high,  black-body-like 
emissivity  of  foam  (Section  5.1).  The  high  foam  emissivity  at  microwave  frequencies  is  explained  with 
the  harmonious  mix  of  contributions  from  all  radiative  processes  taking  place  in  foam.  There  are  two 
necessary,  closely  linked  components  that  ensure  foam  high  emissivity,  namely  high  absorption  losses 
in  a  coupled  foam-seawater  system  and  effective  transmission  of  electromagnetic  radiation  into  and 
through  the  foam.  We  show  that  the  most  important  reason  for  the  high  foam  emissivity  is  the  vertical 
profile  of  the  foam  properties.  Only  in  vertically  structured  foam  layers  does  every  property  and 
process  contribute  toward  delivering  maximum  radiation  to  the  most  effective  absorptive  agent  thus 
producing  maximum  absorption  and  emission. 
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The  conceptual  understanding  of  the  properties  and  processes,  which  determine  the  sea  foam  emissivity 
gives  physical  foundation  and  facilitates  the  implementation  aspects  of  a  foam  emissivity  model. 
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Appendix:  Scattering  Regimes  Mapped  in  m  -  jc  Domain 

For  convenience,  this  appendix  summarizes  information  on  scattering  regimes  relevant  for  this 
study  using  the  work  of  van  de  Hulst  ([37],  Section  10.1).  The  Mie  solution  of  the  Maxwell  equations 
gives  rigorous  mathematical  representation  of  the  scattering  and  absorption  of  EM  radiation  by 
particles.  In  Mie  theory,  absorption,  scattering  and  extinction  efficiency  factors  {Qa,  Qs,  and  Qe,)  are 
formulated  in  terms  of  size  parameter  x  and  refractive  index  m.  Values  of  x  could  be  from  0  to  0°;  m 
varies  from  1  to  0°  for  scatterers  in  vacuum  or  could  be  less  than  1  in  other  surrounding  media. 

Complete  computations  for  arbitrary  combinations  of  x  and  m  values  are  laborious,  thus  limiting 
cases  of  x  and  m  are  usually  investigated  and  approximate  computational  methods  are  offered. 
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Available  approximate  methods  (and/or  theories)  for  various  seattering  oases  are  oonveniently 
presented  in  the  so-oalled  m-x  domain  ([37],  Seotion  10.1;  [52],  p.  408).  Various  seattering  regimes  are 
defined  by  the  values  of  the  seattering  parameters  following  these  inequalities  ([37],  pp.  132-133): 


x«  1 

or 

x»  1 

(A.1) 

m  —  1 «  1 

or 

m  —  1»  1 

(A.2) 

X  (m  -  1)  «  1 

or 

X  (m  -  1)  »  1 

(A.3) 

For  oomplex  m,  oondition  A.2  is  |m  -  1|  «  1  and  it  means  that  both  Re{m}  -  1  and  |Im{m}|  should 
be  «  1  ([37],  p.  179).  Following  these  oonditions,  various  seattering  regimes  are  assigned  to  different 
portions  of  the  m-x  domain.  These  seattering  regimes  (aoronyms  in  Table  A)  are  shown  in  Figure  9, 
whioh  presents  a  simplified  version  of  the  m-x  domain  given  by  van  de  Flulst  ([37],  Figure  20]). 

Table  A,  Approximate  seattering  theories  for  boundary  regions  of  the  m-x  domain 
(modified  from  Table  8  in  van  de  Hulst  ([37],  p.  133). 


Region 

X 

m-1 

x{m  -  1) 

Regime 

RG(1) 

arb 

s 

s 

Rayleight-Gans 

AD  (2) 

1 

s 

arb 

Anomalous  Diffraction 

GO  (3) 

1 

arb 

1 

Geometric  optics 

TR  (4) 

arb 

1 

1 

Total  Reflector 

OR  (5) 

s 

1 

arb 

Optical  Resonance  (Mie) 

RS  (6) 

s 

arb 

s 

Rayleigh  scattering 

The  upper  limits  for  mierowave  frequeneies  in  the  m-x  domain  (Figure  9)  are  as  follows.  For  the 
size  parameter,  x  >  50  is  possible  ([52],  p.  408).  The  largest  value  of  Re{m}  eneountered  in  praetiee 
eould  be  that  of  pure  water,  approaehing  10  asymptotieally  for  frequeneies  below  10  GHz  ([61],  Figure  2). 
As  for  |Im{m}|,  the  high  values  of  the  imaginary  part  of  seawater  permittivity  [62],  Figure  3  (eited  also 
by  Maul  [63],  Figure  5.4),  keep  its  upper  limit  eonsistently  high  above  4  (Table  3)  for  frequeneies 
below  100  GHz. 
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